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ABSTRACT  We  evaluated  fragments  of  the  mitochondrial  COI  and  ND6  genes  to  explore  phylo¬ 
genetic  relationships  among  13  of  the  20  species  of  the  Leucosphyrus  Group  of  Anopheles  (Cel lie) 
(Diptera:  Culicidae),  including  all  four  of  the  currently  recognized  complexes.  Nucleotide  sequence 
data  were  analyzed  using  maximum  parsimony,  maximum  likelihood,  and  Bayesian  methods.  The 
results  revealed  the  monophyly  of  the  Leucosphyrus  Group  and  the  Hackeri  and  Riparis  Subgroups; 
however,  the  Leucosphyrus  Subgroup  and  the  Leucosphyrus  Complex  were  recovered  as  polyphyl- 
etic.  The  monophyly  of  the  Dims  Complex  was  corroborated  by  all  the  analyses  but  with  discordance 
in  the  placement  of  An.  balabacensis  Baisas.  The  maximum  parsimony  strict  consensus  tree  and 
maximum  likelihood  topology  support  the  placement  of  An.  balabacemis  within  the  Dirus  Complex, 
whereas  the  Bayesian  topology  placed  the  species  as  sister  to  the  Hackeri  and  Riparis  clade.  Support 
for  the  split  leading  to  An.  Uitens  Sallum  &  Peyton  and  An.  leueosphyrtts  Donitz  is  not  strong;  however 
in  the  maximum  likelihood  topology  by  using  PHYML,  they  were  recovered  in  a  basal  group  within 
the  Leucosphyrus  Group. 
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The  Leucosphyrus  Group  belongs  to  the  Neomyzo- 
myia  Series  of  Anopheles  ( Cellia )  (Diptera;  Culicidae) 
(Harbach  2004;  Sallum  et  al.  2005a, b)  and  includes  20 
named  species  and  two  geographical  forms  (Peyton 
1989).  Six  species  of  the  Leucosphyrus  Group  are  of 
great  epidemiological  importance  as  highly  compe¬ 
tent  vectors  of  human  malaria  parasites  in  Southeast 
Asia:  Anopheles  balabacensis  Baisas  (White  1983, 
Schultz  1992,  Barcus  et  al.  2002),  Anopheles  latent 
Sallum  &  Peyton  (Zulueta  1956,  White  1983),  Anoph¬ 
eles  leucas])hijrus  Donitz  (Warren  et  al.  1963),  Anojih- 
eles  baimaii  Sallum  &  Peyton  (Rahman  et  al.  1977, 
Rosenberg  and  Maheswary  1982,  Dutta  et  al.  1991, 
Prakash  et  al.  2001),  Anopheles  dints  Peyton  &  Harri¬ 
son  (Eyles  et  al.  1964;  Scanlon  and  Sandhinand  1965; 
Sloof  and  Verdrager  1972;  Ismail  et  al.  1974,  1975; 
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Wilkinson  et  al.  1978;  Deng  et  al.  1982;  Trung  et  al. 
2004),  and  Anopheles  stdawesi  Koesoeniawinangoen 
(Warren  and  Wharton  1963).  Other  species  of  the 
group  are  suspected  to  transmit  simian  malaria  para¬ 
sites  (Warren  and  Wharton  1963,  Coatney  et  al.  1971, 
Tsukamoto  et  al.  1978,  Fooden  1994). 

The  current  classification  of  the  Leucosphyrus 
Group  was  initially  proposed  by  Colless  (1956)  and 
Reid  (1968)  and  later  corroborated  by  Peyton  (1989). 
Subsequently,  Peyton  proposed  the  Elegans,  Leuco¬ 
sphyrus  and  Riparis  Subgroups  based  on  morpholog¬ 
ical  similarities.  The  Leucosphyrus  Group  was  dem¬ 
onstrated  to  be  monophyletic  and  the  earliest 
diverged  lineage  within  the  subgenus  Cellia  (Sallum  et 
al.  2000).  Species  of  the  Leucosphyrus  Group  were 
defined  mainly  based  on  morphology  (for  details,  see 
Sallum  et  al.  2005b),  but  the  12  species  included  in  the 
Leucosphyrus  Subgroup,  plus  An.  tnirans  Sallum  & 
Peyton  (Hackeri  Subgroup),  were  investigated  using 
a  multidisciplinary  approach  that  included  morphol¬ 
ogy  (Peyton  and  Harrison  1979,  Peyton  and  Ranral- 
ingarn  1988),  karyotypes,  polytene  chromosomes,  and 
crossing  studies  (Baimai  et  al.  1984a,b,  1987. 1988a,b,c; 
Bairnai  and  Green  1985;  Sawadipanich  et  al.  1990; 
Poopittayasatapom  and  Baimai  1995).  Consequently, 
to  distinguish  among  the  species  it  is  necessary  to  use 
all  life  stages  (Sallum  et  al.  20Q5a,b),  ultrastructure  of 
the  eggs  (Damrongphol  and  Baimai  1989),  and  alter¬ 
native  idcntiii cation  methods  such  as  those  of  Baimai 
et  al.  (1987, 1988b,c), Sawadipanich  et  al.  (1990),  Wal¬ 
ton  et  al.  (1999),  Huong  et  al.  (2001),  and  Manguin  et 
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al.  (2002).  The  Leucosphyrus  Subgroup  includes  the 
Dims  Complex,  the  I^eucosphyrus  Complex,  the  un- 
iissigned  Anojdieles  baisasi  Collcss,  and  the  geograph¬ 
ical  Con  Son  Form.  The  Dims  Complex  comprises 
seven  species,  An.  dims.  Anopheles  cracens  Sallum  & 
Peyton,  Anopheles scanltmi Sallum  &  Peyton,  An.  baim- 
aii.  Anopheles  nemophilotts  Peyton  &  Ramalingam, 
Anopheles  clegans  (James),  ;uid  Anopheles  lakasagoen- 
sis  Morishita.  The  Leucosphyms  Complex  includes 
An.  laicosfdujrus.  An.  lalens,  Anopheles  introlalns  Col- 
less,  and  An.  balabacensis.  The  Riparis  Subgroup  con¬ 
sists  of  Anopheles  riparis  King  &  Baisas,  Anojdieles 
crisiattis  King  &  Baisas,  Anopheles  macarthuri  Colless, 
and  the  Negros  Form.  Sallum  et  al.  (2005a)  transferred 
An.  elegans  to  the  Dims  Complex  and  thus  renamed 
the  Elegans  Subgroup  as  the  Hacked  Subgroup  to 
reflect  the  change.  Currently,  the  Hackeri  Subgroup 
includes  An.  hackeri  Edwards,  An.  pujutensis  Colless. 
An.  mirans.  An.  Sulawesi,  and  Anopheles  recens  Sallum 
&  Peyton. 

Several  studies  using  genetic  and  molecular  tools 
were  carried  out  to  investigate  species  recognition, 
gene  flow,  and  genetic  population  stmeture  of  mem¬ 
bers  of  the  leucosphyms  Croup  (for  review,  see  Sal¬ 
lum  et  al.  2005b),  but  few  studies  have  addressed 
phylogenetic  relationships  among  members  of  the  Di¬ 
ms  Complex.  Crossing  experiments  (Baimai  et  al. 
1987),  cytologieal  studies  (Baiinai  et  al.  1980,  1988c), 
and  allozymc  analysis  (Green  et  al.  1992)  all  suggested 
a  sister  relationship  between  An.  dims  and  An.  scan- 
luni,  with  An.  bainuiii  being  more  distantly  related.  In 
contrast,  Walton  et  al.  (2000,  2001)  found  that  An. 
dims  and  An.  baiinaii  are  genetically  more  similar  to 
each  other  than  to  An.  dims  or  An.  scanloni.  More 
recently,  Manguin  et  al.  (2002)  observed  that  An. 
scanloni  shares  sequence  characterized  amplified  re¬ 
gions  (SCAR)  fragments  with  An.  dims. 

The  objectives  of  this  study  were  1)  to  test  the 
monophyly  of  the  Leucosphyms  Group;  2)  to  test  the 
monophyly  of  the  Leucosphyms,  Riparis,  and  Hackeri 
Subgroups;  and  3)  to  estimate  phylogenetic  relation¬ 
ships  among  taxa  of  the  leucosphyms  Group.  Two 
geographical  Forms,  Con  Son  and  Negros,  were  not 
included,  nor  were  seven  other  species  for  which  we 
could  get  neither  fresh  specimens  nor  DNA  from  mu¬ 
seum  specimens. 

Materials  and  Methods 

Collection  data  are  in  Table  1 .  In  this  study,  we  used 
fragments  of  the  mitochondrial  cytochrome  c  subunit 
I  (COI  mtDNA)  and  the  NADH  dehydrogenase  sub¬ 
unit  six  genes  (ND6  mtDNA)  derived  from  museum 
specimens.  Leucosphyrus  Group  species  identifica¬ 
tions  were  confirmed  by  either  moiphology  or  poly- 
tene  chromosomes  (for  details,  see  Sallum  et  al. 
2005b).  All  specimens  are  deposited  in  the  Smithso¬ 
nian  Institution,  National  Museum  of  Natural  History 
(NMNH)  collection.  Most  adult  specimens  were  in¬ 
dividually  reared  with  fourth  instar  larval  and  pupal 
exuviae,  and  adult  male  genitalia  kept  as  vouchers. 
When  possible,  we  used  progeny’  brood  specimens 


that  originated  from  individual  wild-caught  adult  fe¬ 
males  subsequently  identified  by  V.  Baimai  by  using 
polytene  chromosomes  (for  details,  sec  Sallum  et  al. 
2005b).  The  remaining  individuals  arc  stored  dry  in 
the  NMNH  where  they  remain  at  ambient  tempera¬ 
ture.  DNA  vouchers  are  stored  at  — 80°C  at  NMNH. 

DNA  Extraction,  Polymerase  Chain  Reaction 
(PCR)  Purification,  and  Sequencing.  Total  DNA  was 
extracted  using  the  DNeasy  tissue  kit  (QIAGEN,  Va¬ 
lencia,  CA)  following  the  manufacturer's  animal  tissue 
extraction  protocol.  DNA  template  was  eluted  in  50  pi 
of  buffer  AE.  Because  the  chance  of  cross-contami¬ 
nation  is  high  when  using  museum  specimens,  we  used 
negative  controls  for  both  DNA  extractions  and  PCRs. 
The  two  primers  used  to  amplify  250  bp  of  the  COI 
gene  were  UEA9.2  (5’-CTA  ACATTTTTTCCTCAA 
CAT  ITT  TTA  GC-3')  and  UEA10.2  (5'-TTA  TTA 
Grr  AAT  AAY  CGT  ART  TCT  C-3' ) ,  both  designed 
for  this  study.  PCR  reactions  were  carried  out  in  a  total 
volume  of  50  pi  by  using  standard  protocols  (Palumbi 
1996).  PCR  amplification  profile  consisted  of  2  min  at 
95°C,  five  cycles  of  1  min  at  94°C,  40  s  at  37°C  and  40  s 
at  72°C,  followed  by  45  cycles  of  40  s  at  94°C,  40  s  at 
48°C,  and  40  s  at  72°C.  PCR  amplification  was  termi¬ 
nated  with  an  extension  of  7  min  at  72°C.  The  two 
primers  used  to  amplify  349  bp  of  the  ND6  gene  were 
ND6.F2  (5'-TTG  GWC  GTA  AWG  GWC  CAT  AAA 
A-3' )  and  ND6.R3  (5'-CAR  GAA  TYT  ATG  TAA  AAA 
CAT  11 1  G-3' ) ,  both  also  designed  for  this  study.  PCR 
reactions  were  performed  under  similar  condition  to 
COI  gene.  The  thermocycling  profile  consisted  of  one 
cycle  of  2  min  at  94°C,  five  cycles  of  1  min  at  94°C,  40  s 
at  37°C  and  40  s  at  72°C,  followed  by  -15  cycles  of  45  s 
at  94°C,  45  s  at  50°C  and  1  min  at  72°C,  with  a  final 
extension  of  7  min  at  72°C.  PCR  products  were  elec- 
trophoresed  in  2%  Tris  borate-EDTA  agarose  gels 
stained  with  ethidium  bromide.  PCR  products  were 
cycle  sequenced  in  both  directions  after  further 
cleanup  by  using  polyethylene  glycol  (PEG)  precip¬ 
itation  (20%  PEG  8000  and  2.5  N  NaCI).  The  cycle 
sequencing  reaction  had  a  total  volume  of  10  pi  and 
included  10  pmol  of  each  primer  and  1  /xl  of  Big  Dye 
terminator  version  3.1.  The  sequencing  reaction  pro¬ 
tocol  consisted  of  one  cycle  of  1  min  at  96°C  followed 
by  30  cycles  of  10  s  at  96°C,  5  s  at  55°C  and  4  min  at 
60°C.  Sequences  were  analyzed  on  an  ABI  Prism  3100- 
Avant  Genetic  Analyzer  (Applied  Biosystems,  Foster 
City,  CA).  The  COI  and  ND6  sequences  were  trans¬ 
lated  into  amino  acids  by  using  the  Drosophila  genetic 
code  implemented  in  MacClade  4.0  (Maddison  and 
Maddison  2000)  and  rcehecked  to  ensure  that  there 
were  no  frame  shifts.  The  sequences  have  been  depos¬ 
ited  in  GenBank  (COI,  accession  nos.  DQ897936- 
DQ897972;  ND6,  accession  nos.  DQ899796  - DQ899832 ) . 

Phylogenetic  Analysis.  Unweighted  parsimony  was 
performed  in  PAUP  (Swofford  2004)  by  using  a  heu¬ 
ristic  search  with  trce-bisection-reconnection  (TBR) 
and  1,000  random-taxon  additions.  Parsimony  boot¬ 
strap  support  values  were  generated  from  1,000  pseu- 
dorcplicates  with  10  random-taxon-addition  replicates 
per  pseudoreplicate.  Parsimony  uninformative  char¬ 
acters  were  excluded  from  all  the  analyses. 


Table  1.  Lorn  lion,  (axon  ID,  «)atr,  rollcrlor*  ami  of  the  »aiuple»  of  the  »perie»  of  the  Leiico»phyni»  Group  included  in  the  study 


Species 

Taxon  ID 

Dale 

Collector 

Geographical  coordinates 

Sex 

Country' 

Location 

An.  lalairl 

IDK43-10 

I8-2.3-K-19S6 

R.  Harbach  and  V.  Baimai 

3°  49'  S  115°  '  E 

M 

Indonesia 

Kalimantan,  Tana  laiut,  SiJaman 

An.  laU-ns\ 

PB53L-100 

I2-XII-19S6 

V.  Daimni 

6°  4’  N  100“  I’  E 

F 

Thailand 

Songkhala,  Sadao.  Pad an g  Besa.  Khao  Rnp  Chang 

An.  tcucasjthynisl 

ID1-002-13 

7-IV-1986 

R.  Harbach  and  V.  Baimai 

1°  2'  S  102“  '  F- 

M 

Indonesia 

Sumatra  Island,  Propinst  Jambi,  Bukit  Bara  (near 

Muarabungo) 

An  It'ucaspht/rusZ 

ID1-007-I I 

7-1V-1986 

R.  Harbach  and  V.  Daimni 

1“  2’  S  102“  '  E 

M 

Indonesia 

Sumatra  Island,  Propinsi  Jambi,  Bukit  Bara  (near 

Muarabungo) 

An.  halahactiwisl 

MI7 

21-X1-1096 

R.  Harbach 

5“  F  S  117“  45'  E 

F 

Malaysia 

Sabah.  Lahad  Datu  District,  Lahad  Datu.  Borneo 

rainforest 

An.  iMilabarvnxisrZ 

Coll.#53 

1989 

M.  Bungs 

3°  51'  S  115“  13'  E 

M 

Indonesia 

south  Kalimantan,  Salaman,  Kintap,  kilometer  IS 

An.  luilah<icenMi\^ 

M47-I6 

23-XM996 

R.  Harbach 

5“  1'  S  117“  45'  E 

F 

Malaysia 

Sabah,  Lahad  Datu  District.  Lahad  Datu,  Borneo 

rainforest 

An  r/ims'3 

09147(29)-3 

4-VU1-19S2 

AFRIMS 

14“  16'  N  101“  54'  E 

M 

Thailand 

Prachinburi,  Ban  Bti  Phrurn 

Ail.  ilimtA 

TH1746(6)-17 

4-X-I9S9 

V.  Baimai 

16“  40'  N  98“40'  E 

M 

Thailand 

Tak.  Mae  Sot,  Ilium  Run 

An.  tlirtitS 

B(12)l 

M 

Thailand 

Tliailand.  Bangkok  colony 

An.  dims (5 

D(15)-ll 

M 

Thailand 

Thailand,  Bangkok  colony 

Ail.  criKtml 

MH00)6(l>-4 

28-1V-1982 

V.  Baimai  and  R.  G.  Andre 

5“  4'  N  102“  56'  E 

M 

Malaysia 

Terengganu,  Kampong  Dura 

An.  crwmiy2 

MH0023(2)-6 

29-I\'-1982 

V.  Baimai  and  R.  G.  Andre 

5“  6'  N  102“  55'  E 

M 

Malaysia 

Terengganu.  Kampong  Tapali 

An.  \amtoni2 

Gass26(5)-21 

1987 

V.  Baimai 

14“  25'  N99“  17'  E 

M 

Thailand 

Kanchanulniri,  Sai  Yok,  Phu  Toei 

An.  scnn/ouM 

9J20-1SO(5)-S 

1- VIM  981 

AFRIMS 

14“  27'  N  99“  S'  E 

M 

Thailand 

Kanchanaburi,  Tha  Kradan,  Ban  Phu  Taka,  Mu  3 

An  ycyih/ouio 

8-31 

1981 

M 

Thailand 

Kanchanuburi  Colony 

An  baimaii'2 

THI690(10)-10 

7-V111-I9S9 

V.  Baimai 

16”  40'  N  98“40'  E 

M 

Thailand 

Mae  Sot.  Ban  Kuriaug,  Tliuni  Rua 

An  baimaift 

08623(2)-109 

26- VI-1982 

AFRIMS 

8“  17'  N  98“  23'  E 

M 

Thailand 

Phangnga,  Khao  Puk  Chuung  (Chong) 

Ail.  haimaiU 

TH.504<2)-M 

23-11-1986 

V.  Baimai 

8“  35'  N  98“  32'  E 

M 

Thailand 

Phangnga,  Ban  Bang  Kaeo 

An.  ItaiinaiiS 

TKK-A 

1987 

M.  M.Thu  and  Myo-Paing 

17”  23'  N  96“  3'  E 

M 

Myanmar 

Pegu  Division.  Taikkyi.  50  miles,  north  of  Yangon 

An.  tutitnaiiG 

636- l-L 

2-X1-1975 

R.  Rosenberg 

24“  l'N  91“  24'  E 

M 

Bangladesh 

Chaklapimgce,  Forest  Beal 

An  elcgitns  1 

F13(6)-27 

11-14- VII1-198I 

H.  Dim! 

14”  19'  N  75”  r  E,  14“  13'  X  75“  1'  E 

M 

India 

Karnataka.  Shimoga,  Kondagalale  and  Shimoga. 

Keladi 

An.  i'/c"om3 

F13(9)-25 

11-14-VII1-1981 

H.  Bhat 

14“  19'  N  75°  7'  E,  14“  13'  N  75”  1 '  E 

M 

India 

Karnataka.  Shimoga.  Kondagalale  and  Shimoga. 

Ketadi 

An.  ncnntfdtUoittl 

0816S-I1 

26-V-I980 

AFRIMS 

8”  36'  N  98”  32'  E 

M 

Thailand 

Phangnga,  Bun  Bung  Ra  Ko 

An.  iu7ito;i/ii/oii.v3B 

08117(I)-100 

16-X1-I979 

AFRIMS 

14“  25'  N  98“  53'  E 

M 

Thailand 

Kanchanaburi.  Huey  Sai  Yok 

An.  iM’inop/ii/oiis-I 

TII49S-IOO 

28-V-1987 

AFRIMS 

8“  35'  N  98°  32'  E 

M 

Thailand 

Phangnga.  Ban  Bung  Kaeo 

Ail.  /afonigoemisl 

F118(6)-8 

28-IV-I980 

J.  C.  Lein/ AFIUMS 

23“  00'  N  120“  00'  E 

M 

China 

Pciynan,  Tungho,  Taiwan  colony,  subcolony 

Bangkok 

An.  ///fcfl-v<7£oefi*ts2 

FI  18(4  >- 16 

28-IV-19S0 

J.  C  Lein/AFRIMS 

23°  00'  N  120”  00'  E 

M 

China 

Pei  yuan,  Tungho.  Taiwan  colony,  subcolony 

Bangkok 

An.  /flAw^oi'n.sis3 

FI 18(3)-l 

28- IV- 1980 

J.  C  Lein/AFRIMS 

23“  00'  N  120“  00'  E 

M 

China 

Peiyuan,  Tungho.  Taiwan  colony,  subcolony 

Bangkok 

Ah.  in  frail*  1 

124-101 /ACC5I0 

7- VIII- 1 975 

E.  L  Peyton  and  Y.-M.  Huang 

6“  50'  N  80“10’  E 

M 

Sri  Lanka 

Western,  Colombo,  l^abugania  Reservoir 

An  infraiiA'3 

689-100 

22-XII-1977 

11“  25'  N  76“  30'  E 

M 

India 

Madras  [Tamil  Nadu],  Nilgiris,  Buliar 

An  .vu/mmi 

0016-16 

23-K-19S5 

J.  (lii 

0“  35'  S  123“  54'  E 

M 

Indonesia 

Toraul,  Bone-Dtimoga  Forest  Reserve 

An.  mmartlmril 

SL47-I17 

22-111-1965 

AFRIMS 

6“  54’  N  100“  15’  E 

M 

Thailand 

Songkhla,  Ton  Nga  Chang  Waterfall 

Aii.  nutatrtliuri'2 

0.8161-24 

25- V- 1980 

AFRIMS 

8°  .35'  N  98“  32’  E 

M 

Thailand 

Phanpiga.  Ban  Bang  Kaeo 

An.  m/u'itiilmrfi 

TH485- 1 00/  Ace  1 269 

26-V-I9S7 

AFRIMS 

8“  38'  N  98“  32'  E 

M 

Thailand 

Phangnga.  Ban  Bang  11a  Ko 

An.  nmcwlhurio 

HKO15-104 

2004 

AFRIMS 

Not  specified 

F 

Thailand 

Thailand 

An.  mncarlhuriB 

CPOOl-lOO 

2004 

AFRIMS 

Not  specified 

F 

Thailand 

Thailand 

5§ 
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Modeltest  version  3.6  (Posada  and  Crandall  1998) 
was  used  to  choose  a  model  using  the  Akaikc  Infor¬ 
mation  Criterion  (AIC).  This  is  similar  to  the  model 
choice  strategy  used  in  Nylander  et  al.  (2004).  Con¬ 
sequently,  maximum  likelihood  topology  was  con¬ 
structed  under  the  HKYIG  model  by  using  PHYML 
version  2.4.4  (Guindon  and  Gascuel 2003). Support  for 
each  elade  generated  for  unpartitioned  data  sets  was 
assessed  by  100  bootstrap  replicates  by  using  PHYML 
version  2.4.4. 

MrBayes  version  3.0B4  (Huelsenbeck  and  Ronquist 
2001 )  was  used  for  Bayesian  phylogenetic  analyses  of 
the  partitioned  (ND6  position  1 ,  ND6  position  2,  ND6 
position  3,  COI  position  1,  and  COI  position  3)  data  set 
without  COI  position  2.  Modeltest  version  3.6  was  used 
to  choose  a  model  for  each  partition  separately.  Mr¬ 
Bayes  version  3.0B4  does  not  implement  all  models 
suggested  by  Modeltest;  therefore  when  a  subset  of 
the  GTR  model  was  suggested  by  Modeltest  then  the 
GTR  model  was  used  in  MrBayes,  with  the  same 
aniong-site  rate  variation  modeling. 

As  outgroups,  we  used  sequences  from  GenBank 
(COI:  An.  aquasalis  AF417697,  An.  albinunuts 
AF4 17695,  An.  gamhiae  1,20934,  and  An.  quadrimacti- 
latusA  NC  000875)  and  (ND6:  An.  at ittasalis  U35260, 
An.  albimanus  U 35259,  An.  g ambiae  L20934,  and  An. 
quadrimactilalusA  NC_000875). 


Results 

Sequences  of 599  bp  (250  bp  for  COI  and  349  bp  for 
ND6)  from  41  individuals  (four  outgroup,  37  ingroup) 
were  obtained  from  the  mitochondrial  COI  and  ND6 
genes  of  13  ingroup  and  four  outgroup  taxa.  Individ¬ 
uals  with  identical  sequence  were  combined  and  re¬ 
named  to  give  32  unique  sequences.  Identical  se¬ 
quences  are  represented  on  Figs.  1  and  2  as  follows:  An. 
leticosphtjrusl.An.  leucosphynts2  is  “hntcosplujrwtljl"; 
An.  dintsi,  An.  dirus5.  An.  r/irrijG,  An.  baimaii 3,  An. 
baimaii  four  is  “dirus4_5_6_baimaii3_4”;  An.  clegans  1, 
An.  etegans 3  is  “t’Zeg«ml_3";  An.  takasagoensisl,  An. 
takasagqensis2  is  "  takasagoensisl  _2"\  and  An. 
inacarthuriS.  An.  macartlniriS,  An.  macartliurib  is 
“niffcflrt/i«n'3_5_6.”  Consequently,  we  analyzed  32  se¬ 
quences  of  13  ingroup  and  four  outgroup  taxa.  The 
number  of  sites,  constant  and  variable  sites,  and  par¬ 
simony  informative  sites  are  listed  in  Table  2.  Because 
COI  codon  position  2  has  only  two  variable  sites,  it  was 
excluded  from  all  analyses. 

The  best  model  for  ML  and  Bayesian  analyses  was 
chosen  with  the  aid  of  Modeltest,  which  suggested  the 
HKYIG  model.  Modeltest  does  not  test  site-specific 
models  for  partitioned  data,  so  to  choose  how  to  best 
model  the  among-site  rate  variation  (ASRV).  the  tree 
used  by  Modeltest  was  reevaluated  by  maximum  like¬ 
lihood  (ML)  using  PAUP  with  the  HKY  model  and 
various  partitioning  and  ASRV  schemes  (Table  2) .  The 
IG  ASRV  on  unpartitioned  data  suggested  by  Mod¬ 
eltest  was  much  better  than  no  ASRV,  and  better  than 
a  site-specific  model  based  on  partitioning  the  data  by 
gene.  However,  a  site-specific  model  based  on  parti¬ 


tioning  the  data  into  five  codon  positions  (no  COI 
position  2)  had  a  better  likelihood  than  the  unpar- 
titioned  IG  model,  and  so  this  site-specific  model 
with  five  partitions  was  used  in  Bayesian  analysis. 
Because  the  site-specific  model  is  not  implemented 
in  PHYML  (and  only  a  simple  version  is  imple¬ 
mented  in  PAUP),  the  IG  model  of  ASRV  was  used 
in  PHYML. 

All  maximum  parsimony  (MP)  (data  not  shown), 
ML  (Fig.  1)  and  Bayesian  (Fig.  2)  trees  based  on 
combined  COI  and  ND6  sequences  show  nearly  iden¬ 
tical  topologies  except  for  a  disagreement  in  the  po¬ 
sition  of  An.  balabacensis,  which  arises  either  as  a 
separate  monophyletic  group  within  a  major  clade 
formed  by  members  of  the  Leucosphyrus  Group  (Fig. 
1)  or  in  the  Hackcri/Riparis  clade  (Fig.  2).  The  ML 
topology  recovered  using  PHYML  version  2.4.4  (Fig. 
I)  corroborates  the  monophyly  of  the  leucosphyrus 
Group  and  recovered  five  subclades.  The  first  sub- 
clade  |  L(u, )  includes  An.  lettcosphynis  and  An.  Intern, 
both  members  of  the  Leucosphyrus  Complex,  the  sec¬ 
ond  subclade  (H)  includes  An.  snlaicesi  and  An. 
miram,  which  belong  to  the  Hackeri  Subgroup,  a  third 
subclade  (R)  includes  An.  macarthuri  of  the  Riparis 
Subgroup,  a  fourth  separate  subclade  [  L(Ll)  |  formed 
by  An.  balabacensis  of  the  Leucosphyrus  Complex  and 
a  fifth  subclade  (D)  leading  to  members  of  the  Dims 
Complex  (An.  dints.  An.  baimaii.  An.  clegans.  An. 
takasagoensis,  An.  craceivt,  An.  scanloni,  and  An. 
nemoplttlous).  An.  dints  and  An.  baimaii  clustered  to¬ 
gether  in  a  clade  that  is  sister  to  An.  elegans.  Mono¬ 
phyly  of  the  L(ll)  is  ambiguous  because  An.  latens  and 
An.  leticosphynts  sequences  clustered  together  in  a 
poorly  supported  clade  (68%  ML  bootstrap  value) 
(Fig.  1 ).  In  all  other  analyses  the  relationship  between 
An.  leucosjdnjnts  and  An.  latens  is  unresolved  (Fig.  2). 

In  contrast,  the  Dims  Complex  was  recovered  as  a 
monophyletic  group  in  all  analyses.  Relationships 
among  its  members  were  not  entirely  resolved  and 
varied  according  to  the  method  used  for  the  analyses. 
An.  cracens,  An.  takasagoensis,  An.  clegans.  An.  baimaii, 
and  An.  dints  were  recovered  monophyletic.  Contrast¬ 
ing,  the  three  sequences  of  An.  scanloni  did  not  cluster 
together  in  any  of  the  analyses,  whereas  those  of  An. 
nemopbilotts  clustered  together  in  a  clade  within  the 
Dims  Complex.  Anojdieles  clegans  was  always  placed 
as  sister  to  (An.  dints  and  An.  baimaii )  lineage, 
whereas  An.  takasagoensis  was  recovered  as  sister  to 
the  (An.  elegans,An.  dints,  and  An.  baimaii)  subclade. 
The  phylogenetic  position  of  An.  cracens  is  not  well 


Table  2.  Description  of  ih«  I^u(uk]iliyrun  Group  anil  ND6  uml 
COI  sefpieuce*  utteil  in  the  analysis 


Partition 

Sites 

Constant 

Variable 

Parsimony 

informative 

NDGposl 

iit> 

80 

36 

25 

NDCpos2 

1!6 

00 

20 

13 

NDtipos] 

117 

45 

72 

51 

COFposl 

84 

74 

10 

6 

COlpos2 

83 

81 

2 

0 

COIpos2 

83 

20 

63 

44 
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Fig.  1.  Maximum  likelihood  bootstrap  topology  using  the  PHYML  version  2.1.1,  with  the  HKYIG  model.  The  site  specific 
(SS)  model  is  not  implemented  in  this  program. 


supported  because  this  species  was  placed  cither  its 
sister  to  (An.  lakasagoensis,  An.  elegans,  An.  dirus,  and 
An.  baimaii)  in  both  ML  and  Bayesian  topologies 
(Figs.  1  and  2)  or  within  a  polytoinv  in  the  Dints 
Complex  plus  the  An.  luilahacensis  lineage,  in  the  MP 
strict  consensus  tree  (data  not  shown). 


The  hypothesis  for  monophyly  of  the  Leucosphyrus 
Croup  is  in  ML  (97%)  and  Bayesian  (1.0)  analyses  and 
moderately  well  supported  in  MP  analysis  (89%). 
Monophyly  of  the  Leucosphyrus  Complex  is  not  sup¬ 
ported  by  any  of  the  analyses:  thus,  the  complex  seems 
to  be  polyphyletic  because  it  includes  An.  balabaccn- 
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Fig.  2.  Bayesian  topology  generated  using  MrBayes  3. OH  t.  The  data  were  divided  into  five  partitions,  based  on 
codon  position  (N136  position  1,  NDfi  position  2.  ND6  position  3.  COI  position  1  and  COI  position  3).  The  model  given  to 
each  partition  is  GTH1;  HKYG:  GTRC:  GTRI  and  G1TKG.  respectively.  All  parameters  were  "unlinked"  in  the  different 
partitions.  A  repeat  of  this  analysis  was  made,  and  the  consensus  tree  differed  only  in  that  c/egnnvl_3  and  clints'i  reversed 
positions. 


sis.  Also,  the  low  bootstrap  values  do  not  support  a 
paraphyletic  hypothesis  for  the  Dims  Complex,  and 
the  phylogenetic  position  of  An.  balnbaccnsis  re¬ 


mains  unresolved  either  as  sister  to  the  Dirus  Com¬ 
plex  or  as  outgroup  of  the  (Hackeri  and  Riparis) 
clade. 
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Discussion 

COI  mitochondrial  DN  A  sequences  have  been  used 
for  studying  genetic  population  structure  of  species  of 
the  An.  i lints  complex.  Consequently,  an  almost  com¬ 
plete  absence  of  mtDNA  differentiation  between  An. 
dirus  and  An.  baiinaii  could  possibly  suggest  either 
mtDN  A  historical  in t regression  between  these  species 
or  a  selective  sweep  that  originated  in  An.  bairnaii 
(Walton  et  al.  2000).  Additionally.  Walton  et  al.  ( 1999) 
identified  an  An.  scanloni-An.  baiinaii  hybrid  among 
field-collected  specimens,  showing  that  there  is  a  po¬ 
tential  for  introgression  between  this  species  pair. 
Jiggins  (2004)  investigated  the  association  between 
mitochondria  and  inale-killer  Wolbachia  in  two  spe¬ 
cies  of  butterflies  of  the  genus  Acraea  and  showed  that 
these  parasites  can  reduce  intraspecific  polymorphism 
and  cause  interspecific  introgression  of  mtDNA.  Hy¬ 
pothetically,  a  cause  but  Wolbachia  infection  has  not 
yet  been  observed  in  Southeast  Asian  Anopheles  (Kit- 
tayapong  et  al.  2000).  In  agreement  with  Walton  et  al. 
(2000),  results  of  the  current  study  found  identical 
NDG  and  COI  sequences  for  both  An.  dims  and  An. 
baimaii,  but  there  is  no  evidence  for  introgression  in 
any  other  species.  Additionally,  there  seems  to  be  very 
low  intraspecific  variation  in  both  genes,  and  thus  we 
found  identical  sequences  for  An.  leucosphyrus,  An. 
elegans.  An.  takasagoensis,  and  An.  inacartliuri,  whereas 
except  for  An.  dims  and  An.  baimaii  interspecific  vari¬ 
ation  seems  to  be  higher. 

In  a  combined  analysis  of  the  COI  and  ND6  gene 
regions,  the  traditionally  recognized  Leucosphyrus 
Group  was  found  to  be  a  strongly  supported  mono- 
phyletic  assemblage  within  the  subgenus  Cellia.  How¬ 
ever,  our  results  revealed  that  the  current  classifica¬ 
tion  of  Leucosphyrus  Subgroup  is  composed  of 
unnatural  assemblages.  In  none  of  the  topologies  re¬ 
covered  using  different  methods  of  phylogenetic  anal¬ 
ysis  were  the  Dirus  and  Leucosphyrus  Complexes  re¬ 
covered  as  sisters.  We  also  provide  evidence  that  the 
Leucosphyrus  Complex  is  not  monophyletic  because 
An.  balabacensis  did  not  cluster  with  the  two  other 
species  of  the  subgroup  included  in  our  study,  An. 
latens  and  An.  le.ticosplujms.  It  is  noteworthy  that  An. 
balabacensis  was  recovered  either  in  the  clade  leading 
to  the  Hackeri  and  Riparis  Subgroups  (Fig.  2)  or  as  a 
separate  lineage  within  the  Leucosphyrus  Group  (Fig. 
1).  The  Dirus  Complex  is  a  monophyletic  lineage. 
Also,  the  Riparis  and  Hackeri  Subgroups  were  recov¬ 
ered  as  sister  groups  (Fig.  2)  or  as  a  polytomy  in  the 
Leucosphyrus  Group.  Relationships  between  An. 
balabacensis  and  the  Riparis  and  the  Hackeri  Sub¬ 
groups  are  not  supported  by  a  morphological  hypoth¬ 
esis.  According  to  Sallum  et  al.  (2005b)  the  morpho¬ 
logical  distinction  between  the  Leucosphyrus  and  the 
Dirus  Complexes  is  problematic  because  some  char¬ 
acters  used  to  define  the  limits  of  each  species  complex 
are  polymorphic.  Generally,  members  of  the  Leuco¬ 
sphyrus  Complex  can  be  distinguished  easily  from 
those  of  the  Dims  Complex  in  having  the  accessory 
sector  pale  (ASP)  wing  spot  present  on  veins  C,  sub¬ 
costa  and  R,  and  by  the  absence  of  pale  scales  at  the 


base  of  hindtarsomere  4.  However,  An.  balabacensis  is 
polymorphic  for  these  characters  and  thus  can  overlap 
with  members  of  both  the  Dirus  Complex  and  the 
Leucosphyrus  Complex.  A  sister  group  relationship 
between  An.  balabacensis  and  members  of  the  Hackeri 
and  Riparis  Subgroups  has  no  morphological  support 
and  An.  balabacensis  can  be  separated  easily  from 
members  of  the  two  subgroups  (Sallum  et  al.  2005b). 
The  placement  of  An.  balabacensis  as  sister  to  the 
Dirus  Complex  is  more  concordant  with  a  morpho¬ 
logical  hypothesis  than  as  sister  to  the  (Riparis  and 
Hackeri)  clade.  Additionally,  Kanda  et  al.  (1983)  com¬ 
pared  seven  populations  of  members  of  the  I-euco- 
sphyrus  Group  and  showed  that  An.  balabacensis  was 
not  distinct  from  An.  dirus.  Similarly,  Yong  ct  al. 
(1983)  used  15  gene-enzyme  systems  to  compare 
genetic  diversity  in  An.  dirus.  An.  cracens,  and  An. 
balabacensis  and  found  that  the  three  taxa  were  mono- 
morphic  for  all  15  loci  tested.  Consequently,  it  is  ob¬ 
vious  that  An.  balabacensis  is  genetically  closely  re¬ 
lated  to  members  of  the  Dirus  Complex.  It  is  also 
possible  that  An.  balabacensis  represents  a  widespread 
species  complex  in  Southeast  Asia. 

Although  our  results  support  the  monophyly  of  the 
Dirus  Complex,  relationships  among  its  members  are 
only  moderately  to  weakly  supported  (Figs.  1  and  2). 
Generally,  within  the  complex,  we  can  recognize  two 
major  groups,  one  group  poorly  supported  group  con¬ 
sisting  of  An.  nemo))hilous  and  An.  scanloni  (data  not 
shown)  and  a  second  group  leading  to  (An.  cracens 
(An.  takasagoensis  (An.  elegans  (An.  baimaii.  An.  di¬ 
ms)  ) ) ).  Sequences  of  An.  scanloni  2,  An.  scan(oni4,  and 
An.  scanloni 5  did  not  cluster  together  in  any  of  the 
analyses  by  using  several  methods.  In  the  Bayesian 
analysis.  An.  .s«m/oni4  clustered  within  the  An. 
rtemophilous  clade,  whereas  An.  scanloni 2  and  An.  scan- 
loni5  formed  a  subclade  that  is  sister  to  An.  nemophi- 
lous  plus  An.  scanloni4  (data  not  shown).  Walton  et  al. 
(2000,  2001)  demonstrated  An.  scanloni  to  be  a  well- 
differentiated  species  and  that  the  high  degree  of 
differentiation  between  northern  and  southern  pop¬ 
ulations  of  An.  scanloni  was  suggestive  of  the  presence 
of  two  incipient  species.  Our  results  also  suggest  that 
there  might  be  at  least  two  subpopulations  within  An. 
scanloni  because  An.  scanloniA  clustered  with  An. 
nemophilous,  whereas  An.  scanloni2  and  An.  scanloni 5 
either  clustered  together  (data  not  shown)  or  were 
recovered  in  a  polytomy  within  the  Dirus  Complex 
(Figs.  1  and  2).  Placement  of  An.  elegans  in  the  clade 
consisting  of  (An.  dims.  An.  baimaii)  might  support 
the  hypothesis  of  Sawadipanichetal.  (1990)  that  there 
arc  cytogenetic  and  crossing  evidence  that  An.  elegans 
is  an  incipient  sibling  species  of  An.  baimaii. 

Results  of  the  current  study  and  other  studies  on 
genetic  relationships  among  members  of  the  Dinis 
Complex  do  not  always  coincide.  For  example.  Man- 
guin  et  al.  (2002)  showed  that  An.  dirus  and  An.  scan¬ 
loni  are  closer  to  each  other  than  cither  is  to  An. 
cracens  or  An.  baimaii  because  they  share  an  888-bp 
SCAR  fragment.  Similarly,  results  oflaboratory  cross¬ 
ing  experiments  and  polytene  and  mitotic  chromo¬ 
some  (Baimai  et  al.  1987)  are  consistent  in  demon- 
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strating  An.  dints  anti  An.  sctmloni  to  be  closely  related, 
and  that  An.  dirus  and  An.  baimaii  were  genetically  the 
most  incompatible  in  comparison  with  the  other  spe¬ 
cies  tested.  Moreover,  Poopittayasataporn  and  Baimai 
(1995)  suggested  that  An.  baimaii  may  be  the  most 
basal  species  within  the  Dims  Complex. 

Results  of  the  current  study  do  not  fully  resolve 
relationships  within  the  Leucosphyrus  Group.  They 
failed  to  show  monophyly  of  the  Leucosphyrus  Sub¬ 
group  and  the  phylogenetic  placement  of  An.  balaba- 
censis.  However,  this  study  corroborates  the  mono¬ 
phyly  of  the  Leucophyrus  Group  and  the  Dirus 
Complex  and  shows  an  initial  indication  of  the  mono¬ 
phyly  of  the  Hacked  and  Riparis  Subgroups.  A  more 
extensive  sampling  of  species  within  the  Leuco¬ 
sphyrus  Group  will  be  critical  to  test  the  monophyly 
of  the  Leucosphyrus,  Hacked  and  Riparis  subgroups 
and  to  establish  the  phylogenetic  position  of  An.  hala- 
bacemis  and  the  two  poorly  known  Con  Son  and  Ne¬ 
gros  Forms.  In  addition,  it  will  provide  a  stronger  basis 
for  future  biogcographical  studies  and  co-evoiution- 
ary  studies  of  Anopheles  species  in  relation  to  simian 
and  human  malaria. 
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